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Abstract 

We investigate general relativistic effects associated with the grav- 
itomagnetic monopole moment of gravitational source through the 
analysis of the motion of test particles and electromagnetic fields 
distribution in the spacetime around nonrotating cylindrical NUT 
source. We consider the circular motion of test particles in NUT 
spacetime, their characteristics and the dependence of effective po- 
tential on the radial coordinate for the different values of NUT pa- 
rameter and orbital momentum of test particles. It is shown that the 
bounds of stability for circular orbits are displaced toward the event 
horizon with the growth of monopole moment of the NUT object. In 
addition, we obtain exact analytical solutions of Maxwell equations 
for magnetized and charged cylindrical NUT stars. 



KEY WORDS: Relativistic stars; gravitomagnetic charge; particle mo- 
tion; electromagnetic fields. 
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1 Introduction 



At present there is no any observational evidence for the existence of grav- 
itomagnetic monopole or so-called NUT charge which may be source of 
gravitomagnetic field in static case. However there are several attempts 
and proposals to detect it through astronomical observations in galactic, 
extragalactic range and in Solar system. For example, the anomalous ac- 
celeration of Pioneer satellites is attempted to explain through the effect 
of gravitational field of magnetic mass on motion of test particles 0E1- 
In recent papers EJ , the theoretical effect of magnetic mass in NUT 
space on the microlensing light curve and the possibility of magnetic mass 
detection using the gravitational microlensing technique have been studied. 

It is not less interesting to study vacuum axial-symmetric solutions of 
Einstein equations describing the gravitational field of nonrotating NUT 
source. Nowadays cylindrically symmetric stationary metrics are analyzed 
due to the fact that they can be used as simple first order approximation 
model for describing some astrophysical objects. Here we will concentrate 
in studying properties of NUT sources related to the cylindrical symmetry. 
The NUT metric components generated by the solution to the Newto- 
nian Laplacian equation corresponding to a Newtonian potential of a rod 
of length 2 A and mass per unit length being equal to 1/2 was presented 
in [H| ■ In the paper |Jj the metric of cylindrical NUT source in the ap- 
proximation of weak field has been obtained and studied. In the paper |K| 
we have shown the metric of Nouri-Zonoz describing spacctime outside a 
line gravitomagnetic field could be obtained from Papapetrou solution of 
vacuum equations of gravitational field. 

Physical interpretation of NUT solution was also investigated in the 
literature. In particular, it was shown in recent paper [5] that the sources 
of the NUT spacetime are two semi-infinite counter-rotating rods of nega- 
tive masses and a finite static rod of positive mass. It was also noted that 
this model is the only possibility for the NUT solution to be stationary 
and possess zero total angular momentum. Although the metric is not, 
from the mathematical point of view, cylindrically symmetric, it is phys- 
ically cylindrically symmetric, which is clear from the gravitoelectric and 
gravitomagnetic fields which are both explicitly cylindrically symmetric 
functions. 

In our preceding paper |l()j the electromagnetic fields of conducting 
shell embedded in the weak gravitational field of cylindrical gravitomag- 
netic monopole were investigated and analytical solutions to the Maxwell 
equations in this spacetime were found. Here we extend them to the space 
time of cylindrical NUT source obtained in P] as exact solution of Einstein 
equations without making any assumption on the strength of gravitational 
field. We are interested to study the motion of test particles and electro- 
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magnetic fields in NUT space with the aim to get tools for studying new 
important general relativistic effects which are associated with nondiagonal 
components of the metric tensor and have no Newtonian analogues due to 
their origin. 

The paper is organized as follows. In section [3 we study the motion 
of test particles in spacetime of cylindrical NUT source. We investigate 
the effective potential and consider the influence of NUT parameter on 
characteristics of particles motion. In section[3]we solve Maxwell equations 
in the exterior metric of cylindrical NUT source around (i) magnetized star 
with the dipolar magnetic field and (ii) electric charge put on cylindrical 
gravitomagnetic monopole. We summarize our results in section 

Throughout, we use the units G = c = 1 and the metric signature 
(—,+,+,+). The round and square brackets on the indices denote sym- 
metrization and antisymmetrization respectively, e.g., Atfiv) = {A^v + A vll ) /2 
and A^] = (A^ - /2. 



2 Particles Motion in a Space of Nonrotating 
Cylindrical NUT Source 

The general stationary axially symmetric solution jSJ[!5| to the vacuum Ein- 
stein equations in the standard Weyl-Papapetrou cylindrical coordinates 
{t, p, ip, z) is given by 

ds 2 = f- 1 [e 2l (dp 2 + dz 2 ) + p 2 d<p 2 ] - f (dt - ujd<p) 2 , (1) 

where the metric coefficients /, 7, uj are functions of p and z only. The 
explicit expressions for /, 7, and to have been obtained by Gautreau and 
Hoffman 6 in the form 

A 2 (x 2 -1) 27 x 2 -l 



(Ax+M) 2 +P ' x 2 -y 
P++P- = P+ -P- 

2A ' V ' 2A 



V ± = VP 2 + (z± A) 2 , A = ^M 2 + I 2 , (2) 

where C is an arbitrary real constant. 

Following to [H] we will take here C = (case C = also corresponds 
to the general NUT metric 

ds 2 = ~a[dt-2l coa9d(/)} 2 +sin 2 9(r 2 +l 2 )d(f) 2 +0-^^ + {r 2 +l 2 )d0 2 , (3) 

which can be received from the quadratic form Q with help of the following 
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transformations 

N 2 ~(M 2 + l 2 ) , „, (N + M) 2 + l 2 
f = — - f _1 e 27 = — - — r = N + M 

1 (N + M) 2 +l 2 ' 1 p+p_ ' + 

^ = , cos (9 = P+ 2A P ~ , w=^0 + -p_), 



v/(r-M) 2 -^ 2 sin (9 , z = (r - M) cos 6 , (4) 



used in paper [S]), parameter a = (r 2 — 2Mr—l 2 ) / (r 2 +l 2 ). Here parameters 
M and I represent gravitoelectric (total mass of the source) and gravitomag- 
netic charges respectively. The parameter I is also called NUT parameter 
or the most commonly known as 'gravitomagnetic monopolc'. From as- 
trophysical point of view it is reasonable to assume that the string has no 
electric charge and consequently in the space-time around the nonspinning 
NUT string the metric tensor has the following nonvanishing components 

.9oo = -/ , 5n = 533 = r 1 ^ 1 , 

902 = fw , 522 = /~V - /w 2 . (5) 

The Euler-Lagrange equations 

d\\d±^)~ dx<*~ ' ( ) 

governing geodesies of the metric Q can be derived from the Lagrangian 

dx^ dx v 

where a dot stands for differentiation with respect to an affine parameter 
A. 

The integrals of motion 

dp t dC 

d Pv dC n t fa\ 

~dX = ^ = '^ = L (8) 

follow from Euler-Lagrange equations 10. Here p a = dC/dx a is canonical 
momenta and the conserved quantities represent, respectively, the total 
energy, orbital angular momentum of the test particle . 

Knowing that g a pp a p^ = —m 2 we can conclude that value of C in J7J) 
equals to — m 2 /2. For spacetime (QJ the Lagrangian Q is 

c = -In 2 + fu><H + Jrvv + i 2 ) + krv - f^w , (9) 
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and the corresponding canonical momenta JSH^J are 



E = f(t - uxp) , L = ftot + (r L p 2 - . (10) 

Our task here is complicated due to the fact that the metric coeffi- 
cients depend not only on radial coordinate, but also on coordinate z (the 
spacetime is axially symmetric). We have to express z through p. In the 
paper j!2j the motion of a free uncharged test particle in the Kerr-Newman- 
Taub-NUT spacetime was considered. Solving the Hamilton- Jacobi equa- 
tion which describes the motion of test particle in the vicinity of non ro- 
tating source endowed with gravitomagnetic charge it was shown that the 
orbit will lie on a cone with 



cos 9 = ± 



21E 



(11) 



As I — > such an orbit degenerates to a single orbit on the equatorial plane. 
Using transformations one can get the expression for z in cylindrical 
coordinates 

21E lL 2 {A 2 +p 2 )-Al 2 A 2 E 2 . . 

z = —\j L 2 -4l 2 E 2 ■ (12) 

It is easy to see that z — corresponds to the equatorial plane 9 — tt/2 
in spherical coordinates. In order to investigate the radial dependence of 
effective potential we restrict ourselves to the study of circular geodesies 
neglecting possible motion of particles along the z axis, hence z = 0. This 
implies that the equation © is 

- ft 2 + 2foj<pt + f-^p 2 + (/-V 2 - = -m 2 ■ (13) 

From equations (Hfl) one could have 

. E(p 2 -cj 2 f 2 )+Lu;f 2 . L-Ecj 

t = 7-2 ! <P = ,_i 2 • (I 4 ) 

JP J P 

Consider a case when particle has a mass m^O. Then it is convenient 
to normalize quantities E, L to the unit of mass of the particle m and 
label them with a bar as E = E/m, L = L/m. The 'tilded' quantities 
are normalized to the total mass of the source M. Then the equation of 
circular motion (|13|) can be rewritten in the form 

r 1 -* (s§)' -"<**>• "ta*>-n-T (is) 

The elliptic integrals resulting from the solution of the equations of 
motion are not particularly informative, but one can get a general picture 
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of the orbits by considering behaviour of an effective potential. Define an 
effective potential V e ff as that value of E such that p = at radius p: 



1 V e 2 ff f(L/M-V eff X)* 



f r 
The solution for effective potential 



= . (16) 



-ujf(L/M) ± y/- f(puP-f ?)fP + p(L/M)*(P 

p 2 — f 2 UJ 2 



Ve ff = tu _ nza 



should have a sense of energy and a root which corresponds to the asymp- 
totic condition linip^oo V e ff — 1 has to be chosen. 

Figure 1 illustrates the dependence of the effective potential on radial 
coordinate p/M for a source endowed with mass M and gravitomagnetic 
mass I. 




Figure 1 : The radial dependence of the effective potential for particles with 
nonzero rest mass for different values of angular momentum L. The left 
hand side figure is responsible for the case when the NUT parameter / = 0. 
For the right hand side figure the NUT parameter I — 0.5. Maxima in 
the effective potential indicate unstable circular orbits and minima stable 
circular orbits. Curves for particles with equal angular momentum and 
different gravitomagnetic charges have more monotonous behavior with 
the increase of the value of NUT parameter. 

To simplify further calculations we can make the assumptions concern- 
ing the value of I. We select second order approximation in gravitomagnetic 
monopole moment which can be justified by the fact that at the moment 
there is no any astrophysical evidence for strong gravitomagnetic mass. 
For example, in the recent papers E] the magnetic mass detection us- 
ing the gravitational microlensing technique has been explored and it has 
been evaluated that the minimum observable gravitomagnetic mass to be 
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about lAmeters. The following assumption is made concerning energy E 
and angular momentum L of the particle. As the quantity I 2 is small then 
for the nonrelativistic particles with small values of energy and angular 
momentum one can neglect terms being proportional to E I 2 and L I 2 and 
get 

21E 



. _, y/T+j? . (18) 

The conditions for circular motion of test particle at radius p = po are 

V(p ) = , V'\ P=P0 = , (19) 

where prime denotes the derivative on radial coordinate. Furthermore, the 
circular motion is stable if 

V"\ P = po < . (20) 

From the equation (|19f) after some lengthly algebra one can obtain the 
following expression for the energy of the particle 



E 2 



2f?(X + f?) 


2? {X 2 - 




-p 2 ))+Xp 2 (X + 


P 2 ) 


X 3 


2(H 


-?)(?- 


h^(X 


-2) -4X)+T 2 Y 





(21) 



where X = 1 + yjl + p 2 and Y = 3p 4 (X - 2) + 10X + Sp^X. The circular 

— 2 . 

motion is possible provided that the expression for is nonncgative |13| . 
Circular orbits exist till a limiting case corresponding to the photon orbit 
(E = E/m — ► oo). When NUT parameter is equal to zero this limit is equal 
to p 2 — fPM 2 = 3M 2 . This corresponds to the value r = 3M in spherical 
symmetric Schwarzschild spacetime. The influence of NUT parameter is 
such that the limit of existence of circular orbits decreases. For example, 
for I = 0.5 we have p 2 = 2.49. 

2 

The orbit is maximally bound when E =1. Equating the expres- 
sion l|21l) to unity one can obtain the equation governing the bounds of 
maximally binding orbits. In the field of a cylindrical source endowed with 
NUT parameter, taken as I — 0.1, radius of maximally binding orbit is 
p 2 1=3 7.98. Under the influence of NUT parameter thejralue of maximal 
radius of bound orbits decreases. For comparison: when I = radius of the 
orbit is p 2 = 8 which corresponds to r = AM in Schwarzschild spacetime. 

The orbit situated at an inflection point of the effective radial potential, 
that is, with 

V"\ P=P0 = (22) 

is the least tightly bound and describes the bounds of stability region. The 
solution of this equation for I = is p 2 — 24. This corresponds to r = 6M 
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for spherical symmetric Schwarzschild case. The orbits are stable for p > 
2\/6 and unstable for p < 2^/6. For I = 0.7 we have for limit of stability 
p 2 = 34.41. Therefore, the radius of critical orbit separating stable orbits 
from unstable for nonzero gravitomagnetic monopole increases. 

The influence of the NUT parameter on the motion of a test particle 
can be seen from the following table 



Gravitomagnetic influence on motion 


I 


7? 

r circular 


7? 

t max. bound 


Pstable 





3 


8 


24 


0.05 


2.99 


7.99 


24.02 


0.1 


2.98 


7.98 


24.08 


0.5 


2.49 


7.5 


27.21 


0.7 


1.97 


6.94 


34.41 



where pl ircular defines the radius of last circular orbit, p^ax.bound defines 
the radius of the last bound orbit and p^ tafc ; e defines the radius of first 
stable orbit. 

Hereafter for simplicity we will omit widetilded and overlined labels 
but keep in mind that all quantities are normalized to the mass of the 
source, energy and angular momentum are also normalized to the mass 
of the particle. The four-velocity of the proper observer has the following 
nonvanishing components 

u a = -^jl ,0 ,0 ,o| , u Q = v/jj-1 ,0,c;,oj . (23) 

Components of absolute acceleration w a = u a -pvP of the proper ob- 
server in the gravitational field are given by 

Wa= -yjo ,/, P ,0 . (24) 

The nonvanishing components of the relativistic rate of rotation Ap a = 

U[a,0\ + U[pw a ] are 



A 



k/Jw,,,, A zv = l^]u> <z , (25) 



For the case of zero gravitomagnetic charge all components of relativistic 
rate of rotation are equal to zero. It means that NUT parameter influences 
on the observer as a force which drags the observer and makes him rotate 
around the source. 
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The geodesic equations S 2 x a /Sa 2 = in spacetime (QJ take the form 

+ [2fupp - 2p 2 ojp f - f(p 2 + f 2 Lo 2 )p u ] = , (26) 

P + (7,p - ^)(P 2 - i 2 ) + ^ [ffj 2 + (-2p/ + £„)^] 
+2(7,, - " e~ 2 ->f(Lof, p + fw, p )<pi = , (27) 

^ + (7,, - ^)(i 2 - P 2 ) + ^ [&^ 2 + ffj 2 ] 

-e- 2 if(Lul z + fu, lZ )<pt + 2( 7 ,p - tj)pz = , (28) 



if , ± 
p 2 fp 2 



<P + -Ar/^ + -rn (2/pp - pV,pP + /,**) - f ^) = , (29) 



where 8 denotes the covariant derivative and (if = f^ p p + f. z z, fi^ = co tP p + 
uj, z z, £ p = (p 2 + fuo 2 )f. p + 2/ 8 ww, p , £ z = (p 2 + f 2 u; 2 )f, z + 2f*uw, z . 

For the circular geodesies p — z — only radial (|27|l and z compo- 
nent (|29|l of the geodesic equations are not vanishing 

fj 2 + (~ + r 2 i P W - 2(w/ )P + = , (30) 



F 2 t^ 2 + fj 2 - 2(w/,, + fu,,)<pt = . (31) 

The angular velocity of the test particle is given by w = dtp/dt = 
dtp/da (dt/dcr)^ 1 — ip/i which becomes, using (|3U|1 and (|31|l as 



K± v / K 2 -(f, p + f, z )x 
w = j , (32) 

X 

where 

« = /[/(o;,p + w I ,)+w(/,p + /,«)] (33) 

and 

X = (/,, + U(P 2 + /V) - VP + p + u, s ) . (34) 

Equation (|32|l dehnes the angular velocity of a zero angular momentum 
particle. Due to the effect of dragging of inertial frames the particles are 
dragged by the influence of NUT parameter. This effect weakens with dis- 
tance as 1/p 3 and it makes the gravitomagnetic charge of source measurable 
in principle. 
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3 Electromagnetic Fields of 

Nonrotating Cylindrical NUT Source 

The general form of the first pair of general relativistic Maxwell equations 
is given by 

-Fa/3,7 + ^7«,/3 + ^0y,a = ) (35) 

where F a p is the electromagnetic field tensor expressing the strict connec- 
tion between the electric and magnetic four- vector fields E a , B a . For an 
observer with four- velocity u" obs j , the covariant components of the electro- 
magnetic tensor are given by 

F a p = 2u( obs) [ a E l3 } + ■q aj3l suJ obs - ) B s , (36) 
where ?7 Q /3 7 5 is the pseudo-tensorial expression for the Levi-Civita symbol 

^q/3 7< 5 _ - ^_ e aj3lS , rjap^s = ^/~^ge a p^5 , (37) 

with g = det\g a p\ = —e 4l p 2 f~ 2 for the metric 

A useful class of observers is represented by the "zero angular momen- 
tum observers" or ZAMOs [T2|. These are observers that are locally sta- 
tionary (i.e. at fixed values of r and 9) but who are "dragged" into rotation 
with respect to a reference frame fixed with respect to distant observers. 
They have four-velocity components given by 

K)zamo = b-' (l,0,-^,o) ; «) ZAMO = b(- 1,0,0,0^ , (38) 

where a = p 2 — f 2 oJ 2 , b = p^/JJa. 

The general form of the second pair of Maxwell equations is given by 

F af3 . fi = 4irJ a , (39) 

where J a is the four-current. 

Maxwell equations assume a familiar flat-spacetime form when pro- 
jected onto a locally orthonormal tetrad. In principle such tetrad is ar- 
bitrary, but in the case of a relativistic rotating metric source a "natu- 
ral" choice is offered by the tetrad carried by the ZAMOs. Using i|3~5|) we 
find that the components of the tetrad {e^} = (eg, ep, e^, ei) carried by a 
ZAMO observer are 

el=b- l (l^- f -^,Q\ , (40) 
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e^y-,0,0,1,0 



3? = V/e-^O, 0,0, . 



(41) 
(42) 
(43) 



We can now rewrite the Maxwell equations (|35|) and Ij39(l in the ZAMO 
reference frame by contracting them with Q40 JI — l|43 |l . After some lengthy but 
straightforward algebra, we obtain Maxwell equations in the more useful 
form 



,27 



-B q 



-B 1 



BP) =0, 



(44) 



dB'p f 



(45) 



dt e 2 ^ 



B 



e 1 p , e 7 wf - 
—=E J-B p 



dB £ _ e~tf 
dt Ja 



,p o, 



and 



(47) 



e "'^r E ') + T E ^ + i el f Ei ) =47r 7^^ / ' (48) 



£V«^! +4 ^JP (49) 



-.(cufEP-pB" 



-.(tofE'+pBP) 



p 27 art „2 7 

4^4=^-/0,^), 
/ dt xffa 



(50) 
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(wfE% + pBf 4 



~f dt 



v7 



(51) 



We will look for stationary and axially symmetric solutions of the 
Maxwell equations, taking into account that in the vacuum region around 
the source all components of electric current are equal to zero. We could 
study the realistic configuration of magnetic fields, for example, the dipolar 
magnetic field. Due to the possible smallness of the gravitomagnetic mass 
of the star estimated from some astrophysical observations we may perform 
calculations to the first order in NUT parameter. 



(U 0.8 

s 

«x 
ffl 



a 

■Q. 

ffl 



0.6 



z/M=3. 2 
z/M=2. 2 
e/H=l . 2 
z/M=0. 2 



p/M 



Figure 2: Dependence of general relativistic modification factor 
BgrI ^Newt °f magnetic field on the radial coordinate p/M normalized 
in units of stellar mass. Near to the NUT source the magnetic field will be 
amplified, then in some intermediate region it will be weakened and in the 
asymptotically far zone the influence of NUT parameter is negligible, the 
behavior of field is Newtonian and relation tends to unity. The influence 
of the NUT parameter is more strong near to the source of the z axis. 

According to the transformations Q solution for vector potential in 
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spherical coordinates (see, for example, ^3 1TB] 
3 n sin 2 9 



2M 



( r \ 2 , . 2M\ 1 / r 



can be rewritten in cylindrical ones as 

3M P 2 



8M 3 / 



In/ 



2M(2M + N) 
(M + N) 2 



(52) 



(53) 



where \i is the dipolar magnetic moment. 

This yields the following exact solutions for the nonvanishing compo- 
nents of magnetic field 



B p 



3(j,pe 



:{2MfN tZ {N + 3M) 



&M 3 f(M + N) 3 
+ f. z (M + N) [(M + N) 2 In / + 3M 2 - N 2 ] | 



(54) 



B z = 



3/ie 



8M*f { M + Nj{ 2MfpN ^ N + 3M ^ 
(M + N) 2M(N + 2M)(pf. p - 2f) 

(Af + A0 2 [p/, P (ln/-l)-2/ln/]]} . 



(55) 




Figure 3: The radial dependence of electric field E? for different values of 
the gravitomagnetic monopole I. The effect of the NUT parameter on the 
electric field is becoming important near to the source of the z axis. 

The dependence of relation Bq R / B p Newt on parameter p/M is shown 
on Figure El 
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In the limit of flat spacetime when p = r sin 6, z = r cos 6 and limjwwo = 
the expressions l)54 [l -l)55 [l reduce to the Newtonian ordinary expressions 
for magnetic field of a magnetic dipole |17| 



B" 



M(2^ 2 - P 2 ) 

( 2 2 + p 2)5/2 



B z = 



(z 2 + p 2 ) 5 / 2 ' 



3/izp 



(56) 



Consider now as a toy model the charged NUT star with the monopolar 
electric field configuration 



If Q is the electric charge per unit length of the line tube then the solution 
for radial electric field admitted by Maxwell equations is 



Figure shows that for small values of z and p (near to the source) the 
influence of NUT parameter is noticeable. 

4 Conclusion 

We have investigated the general relativistic effects associated with the 
possible existence of the gravitomagnetic monopole moment of cylindrical 
NUT source and in particular studied the influence of the gravitomagnetic 
monopole moment on the motion of the test particles along circular orbits. 
These effects could be divided into two parts. First effect is related to the 
fact that due to the influence of the NUT parameter the particle's circular 
orbits will be not kept in the equatorial plane. Our aim was mainly related 
to the investigation of the second effect connected with the influence of 
the NUT parameter on the circular orbits of the particles. It is shown 
that under the influence of the NUT parameter the particle's orbits are 
becoming less stable which can be seen from the table presented in the 
paper. 

Analytical general relativistic expressions for the dipolar magnetic field 
around NUT star are presented. It is shown that in the linear approxima- 
tion in NUT parameter there is only effect of mass of the star on the value 
of stationary magnetic field. In the case of charged star there is strong 
influence of NUT parameter on the value of monopolar electric field. Our 
results in principle could be combined with analysis of astrophysical data 
on electromagnetic fields of compact objects in order to get upper limits 
on the value of NUT parameter. 



E" = E"{p,z)^Q , 
E~ z =E* = 0. 



(57) 
(58) 



EP = Q 



(59) 
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